Loss of 7r-Junction Behaviour in an Interacting Impurity Josephson Junction 
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Using a generalization of the non-crossing approximation which incorporates Andreev reflection, 
we study the properties of an infinite-U Anderson impurity coupled to two superconducting leads. 
In the regime where A and Tk are comparable, we find that the position of the sub-gap resonance 
in the impurity spectral function develops a strong anomalous phase dependence- its energy is a 
minimum when the phase difference between the superconductors is equal to it. Calculating the 
Josephson current through the impurity, we find that 7r-junction behaviour is lost as the position of 
the bound-state moves above the Fermi energy. 



What is the Josephson current through an interacting 
Anderson impurity? As this phenomenon involves the 
coherent transport of pairs of electrons through the im- 
purity, one would expect that a strong on-site repulsion 
would greatly diminish the effect. Though this intuitive 
expectation is partially correct, there is a variety of non- 
intuitive behaviour also associated with this system. 

If the on-site repulsion is large and the impurity is 
singly occupied, a lowest-order perturbative calculation 
in the impurity-superconductor coupling reveals that the 
sign of the Josephson coupling can become negative, 
meaning that the corresponding ground state of the sys- 
tem has a phasp|-|diffcrcncc of S = it between the two 
superconductorsEn. An appealing explanation for this 
behaviour was provided by Spivak and Kivelsonn, who 
showed that in the limit of single occupancy, it is im- 
possible to transport a pair across the impurity while 
preserving its spin ordering, leading to an extra factor of 
(—1). More generally, 7r-junction behaviour is epmected 
whenever spin-flip tunneling processes dominatem. Re- 
cent work using a Hartree-Fock type procedure indicates 
this behaviouc can survive beyond lowest-order pertur- 
bation theoryU. 

Alternatively, in the limit of large impurity- 
superconductor coupling, the physics of the Kondo effect 
will also become significant. A resulting enhancement 
of the Josephson current through the impurity has been 
predicted in the regime where the superconducting-gap A 
is much smaller than the Kondo temperature T^-El. Sig- 
nificantly, no 7r-junction behaviour is expected here, as 
the spin of the impurity is completely screened by the 
Kondo effect. 

Given these prior results, it is natural to ask how large 
the ratio Tr-/A must be to see the destruction of ir- 
junction behaviour, and what the properties of the junc- 
tion are in this cross-over region; these questions are 
the motivation for the current paper. We examine the 
regime where A is comparable in magnitude to Tk, a 
regime in which the effects of both Kondo physics and 
superconductivity must be treated on an equal footing. 
Note this parameter range is also of interest as it is 
more appropriate to mesoscopic quantum dot systems, 
which would be one possible experimental realization of 



the model. We find, somewhat surprisingly, that the 
properties of our system can be understood in terms of 
another well-studied feature of magnetic impurities cou- 
pled to superconductors- the sub-gap bound state. Typ- 
ically, rthe position of this state is a function of the ratio 
Tk/Auu. We find now that this resonance also devel- 
ops a pronounced dependence on the phase difference- as 
6 increases from zero, the energy of the bound state de- 
creases. This phase dependence is anomalous in the sense 
that it indicates an energetic preference away from 5 = 0. 
Interpreting this sub-gap state as a current carrying An- 
dreev bound stataL3, we are lead to the conclusion that 
the system is a 7r-junction if the sub-gap bound state is 
located below the Fermi energy Ep. This relation is con- 
firmed by making an explicit calculation of the Josephson 
current through the impurity. 

Formalism. -The model we study consists of an infinite- 
ly Anderson impurity coupled to two superconducting 
leads, each having a different phase. Using a slave-boson 
representation, we have: 

H = H + e d ]T flf a + W J2 (ite&Va + ^ c -), (!) 
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operators here create band electrons, with a 
denoting spin and a = L, R labelling the left and right 
superconducting leads. A Q = |A| exp(i</> a ) represents the 
pair potential in lead a, with the phase difference 6 be- 
ing defined as 4>r — 4>l- The Anderson impurity has bare 
energy Ed, and is represented in the usual manner us- 
ing auxiliary fermion (/) and boson (b) operators; the 
U = oo constraint of single occupancy takes the form 

The NS-NCA.-We calculate the impurity spectral 
function (also called the impurity density of states) for 
our system using the NS-NCALj, an extension qtthe self- 
consistent non-crossing approximation (NCA)li3 to su- 
perconducting systems. The NCA amounts to an in- 
finite order re-summation of perturbation theory, and 
has been shown to be quantitatively reliable down to 



temperatures below Tfs-11313. The modification we em- 
ploy seli-consistently includes multiple-Andreev reflec- 
tion processes in a manner which is formally exact to 
order 1/iV, where N — 2 is the spin degeneracy of the 
impurity. As the success of the normal NC A is attributed 
to the fact that it too is exact to order 1/N (in the ab- 
sence of superconductivity), the NS-NCA employed here 
can be viewed as a natural extension to systems with su- 
perconductivity. The graphs determining the /-fermion 
and b slave boson propagators are given in Fig. 1. 



f-Fermicm Propagator: 



./O. 




b-Slave Boson Propagator 



FIG. 1. Diagrammatic representation of the NS-NCA. 
Dashed lines are f-fermions, wavy lines are slave bosons, solid 
lines are lead electrons. Double lines indicate a fully dressed 
propagator. Anomalous propagators indicate Andreev reflec- 
tion. 

Note that the anomalous graphs appearing in Fig. dfl) 
always involve two Andreev reflections. As these Andreev 
reflection events pick up the phase of the superconductor 
in which they occur, and each Andreev reflection is free 
to occur in either of the two superconductors, the phase 
difference S naturally enters the impurity self energies 
through an interference term. A previous study using the 
NC A to study the regime A > > Tr-EJ neglected these di- 
agrams, and thus omitted the only phase-dependent con- 
tribution to the impurity Green function which survives 
in the limit of a strong on-site repulsion. 

The Dyson equations pictured diagrammatically in 
Fig. 1 lead to a set of coupled integral equations for 
the /-fermion and slave boson propagators. Letting 
F(uj) = [w - e d - EH)" 1 and B(u>) = (w - n(w))- 1 
represent the /-fermion and slave boson retarded propa- 
gators respectively, the equations read: 



E(w) = - /de(p(e)B(w -e)f(-e)- (3) 



- / de'T(e)T(e')B{u) + e)F{w + e + e')B{uj + e') 



n( w ) 



2r 



de(p(e)F(Lu + e)f(e)+ (4) 



de'T{e)T(e')F{uj + e)B{w + e + e')F(u> + e') 



In these equations, p(s) is the electronic density of 
states, T = TrW 2 p(0) the bare tunneling rate, / the Fermi 



distribution function, and r(e) is an effective electron- 
hole coherence parameter defined by: 



T H =^2 u *k, a v k ,J(\ UJ \ - e n )f(u) 



(5) 



k,a 



where Uk t a and Vu,a are the usual BCS coherence factors. 
Results- We use a Gaussian density states for the band 
electrons having half- width D. For the results shown 
here, we choose model parameters Ed = —0.671? and Tj_^p 
0.I5-D, which yields an approximate Tk = 0.0005Tt3. 
We have solved numerically via iteration the NS-NCA 
equations in equilibrium for various temperatures, gap 
sizes and phase differences. Within the NCA, the im- 
purity spectral function Ad(to) can be directly related to 
the /-fermion and slave boson spectral functions: 

A da {uj) = J ds[e-^ s + e-^-^]A fa (e)A b (e - w) (6) 

where the auxiliary particle spectral functions are defined 
by A f = -ilmF, A b = -~ImB. Note that the equal- 
ity in equation (J6|) reflects a neglect of vertex corrections 
which is consistent with the large- N nature of the NCA. 
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FIG. 2. Position of sub-gap bound state as a function of 
A/T K for a fixed T/A = 0.5, T K = 0.0005D. The vertical 
dashed line indicates the approximate value of A for which the 
bound state resonance crosses the Fermi energy. Inset: plots 
of the impurity density of states Ad(u>) for A/Tk — 0.66 at 
various phase differences. From right to left, we have 8 = 
(solid curve), 5 = 7r/4, 5 — it/2, 5 — -k (dash-dot curve). 

The solid curve in Fig. @ (with points indicated by cir- 
cles) shows the position y (in units of A) of the sub-gap 
resonance in the impurity density of states as a function 
of A/Tr- at 5 = and at a constant ratio T/A = 0.5; this 
resonance corresponds to the sub-gap bound state. As 
has been observed in earlier studies, we find that as A/Tk 
decreases, the position of the bound state increases, with 
the transition across the Fermi energy occurring when A 



is on the order of Tk- Note that as we are not using a 
model with particle-hole symmetry, we only have a single 
sub-gap state. 

rA,cogent interpretation of this behaviour was provided 
iriffl, where the position of the sub-gap bound state was 
suggested to result from the competition between two 
possible ground states. For Tk > A, the ground state of 
the system is the Kondo singlet, and the excited sub-gap 
state is a spin doublet. Recall that in the Kondo singlet 
ground state of the m&mte-U Anderson model, there is a 
small probability of finding the impurity empty, despite 
the fact that the bare site energy Ed is well below Ep; 
this is not true for the doublet state, where the impurity 
is expected to be singly occupied. Thus, the fact that the 
bound state is located above Ep in this limit (i.e. y > 0) 
reflects the fact that one must add a particle at the im- 
purity to the singlet ground state to reach the doublet 
state. The magnitude of y indicates the energy splitting 
between these two states: 
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Edoub — Ek 



>0 



(7) 



where Edoub is the energy of the doublet state and Ek is 
that of the Kondo singlet state. Similarly, for A > Tk, 
the ground state is the fully-paired doublet favoured by 
the superconductivity; the bound state is located be- 
low Ep as one must remove a particle from the doublet 
ground state to reach the excited singlet state. An iden- 
tical expression holds for y: 
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doub 



<0 (8) 



The crossing of Ep by the sub-gap state is thus seen to 
indicate a substantial change in the nature of the ground 
state. 

While the behaviour of the bound state at S = has 
received much attention, the effects of having a phase dif- 
ference has not, to our knowledge, been previously stud- 
ied. The triangles in Fig. indicate the bound state po- 
sitions y for the same values of T&-/A as the solid curve, 
but now with a phase difference of 8 — ir between the 
superconductors. We find that for all values tested, y 
decreases as we increase 5 to ir- the sub-gap bound state 
has an anomalous phase dependence. This behaviour is 
seen more explicitly in the inset of Fig. g, where we plot 
the sub-gap resonance for various values of 6. 

The phase dependence of the sub-gap state can also 
be interpreted in terms of a competition between the 
Kondo singlet and the fully paired doublet states. The 
Kondo state is again expected to be the ground state for 
Tk > A, while the the doublet state will be the ground 
state in the opposite regime. Significant now however is 
the fact that the phase dependence of these two states is 
quite different. In the doublet state, the impurity can be 
viewed as a local moment and we thus expect a negative 
Josephson coupling and that Edoub will be a minimum at 
S = 7r. In the Kondo singlet state, the impurity spin has 



been screened and thus we do not expect any anomalous 
Josephson coupling- Ek will be a maximum at 6 = tt. 
With these associations, it follows from Eqs. (Q) and (||) 
that y is minimized at 5 = n, meaning that the sub-gap 
state has an anomalous phase dependence. Note that for 
small values of A/Tk, it is possible to make the sub-gap 
state cross Ep and thus change markedly the nature of 
the system ground state by only changing 5. 

At this point, it is natural to attempt to make a con- 
nection to the Josephson effect. Recall that for a non- 
interacting system, it is possible to discuss the Joseph- 
son current as being at least partially carried by Andreev 
bound states existing in the weak link between the two 
superconductorslMl These states have a phase-dependent 
energy, and thus contribute to the current through the 
relation / = j^dF/d5, where F is the free energy of the 
junction. For the simple case where a non-interacting 
impurity couples the two superconductors, there are al- 
ways two Andreev bound states with energies ±|e(<5)|; as 
the lower-energy state — \s(S)\ has a "normal" phase de- 
pendence (dE/dS > 0), the Josephson current is always 
positive. 

In the present case, matters are quite different due to 
the strong on-site interaction. We have only one appar- 
ent Andreev bound state, and its phase-dependence is 
anomalous (dE/d5 < 0). We would thus naively expect 
that the Josephson current would both become negative 
and enhanced in magnitude as the bound state crosses 
below Ep and becomes "occupied" . Of course, this pic- 
ture is oversimplified, as there are other possible contri- 
butions to the Josephson current- due to interactions, 
the impurity contribution to the free energy F is not just 
a simple function of the impurity density of states. 

The conjecture made above can be tested by making an 
explicit calculation of the DC Josephson current throuffh, 
the impurity. Using a technique similar to that used int£l 
to calculate the normal current through an interacting 
impurity, we arrive at_the following exact expression for 
the Josephson currentEH: 



4e . 8 W 2 
1 j = — I sin - ) > i. 



duf(u}) x 



(9) 



— [9 k M G d (<«>,*) 



R,12 



Here, g k ' (u>) indicates an anomalous BCS propagator, 
and G d (u>,5) is the anomalous impurity propagator; 
without loss of generality, we have chosen (f>i + <pR = 0. 
We calculate the anomalous impurity propagator 
within—the NCA using the lowest-order contributing 
grapblldld, without making any further approximation. 
To calculate an approximate value of the temperature 
Josephson current, we use Eq. with the Fermi function 
taken at T = 0, but with all spectral functions calcu- 
lated at T = 0.5 A (the NCA is unreliable at T = 013). 
This is a reasonable procedure, as we expect no qualita- 
tive changes in the spectral functions as the temperature 
is further lowered; the sub-gap resonances will only be 



sharpened. Retaining a finite T Fermi function would 
smear the contribution of the sub-gap resonances, and 
thus obscure the behaviour we are interested in. 




FIG. 3. Approximate T — phase-derivative of the 
Josephson current (dlj/dS) at 5 = for various values of 
A/Tk- Note the transition from negative to positive cur- 
rent (indicated by the vertical dashed line) coincides with the 
sub-gap bound state crossing the Fermi energy (see Fig. El). 

Plotted in Figure is the approximate temperature 
phase derivative of the Josephson current (dlj/d5) thus 
obtained as a function of A/Tk, taken at 6 = 0. Com- 
parison to Fig. indicates that 7r-junction behaviour is 
indeed lost when the bound state crosses above the Fermi 
energy, as expected from a non-interacting picture. Sim- 
ilar results hold at other values of S- (dlj/dS) changes 
sign at roughly the same value of A/Tk at which the 
bound state crosses the Fermi energy. It would thus ap- 
pear that the sub-gap bound state is the entity responsi- 
ble for carrying the reversed-sign Josephson current as- 
sociated with 7r-junction behaviour. 

Note that in some respects, our conclusion here resem- 
bles what was found in a study of the Josephson curreat 
through a non-interacting magnetic insulating layered; 
there also, the transition to 7r-junction behaviour is ac- 
companied by an anomalous Andreev bound state cross- 
ing below the Fermi energy. We stress again that the 
present situation is quite different due to the strong on- 
site repulsion- there is no a priori reason guaranteeing 
that the general non-interacting theory relating Ij to An- 
dreev bound states should be applicable to the interact- 
ing system studied here. We find only a single sub-gap 
state in the impurity density of states, whereas in the 
non-interacting case bound states always occur in pairs, 
with one member of each having dE/d<p < 0. 

Also of interest in the present system is the behaviour 
of the superconducting phase of the impurity itself. We 
find that this phase undergoes a shift by 7r as the bound 
state passes through the Fermi energy. When the bound 
state is above Ep, the impurity simply has the average 
phase of the two superconducting leads- at low energies, 
the signs of both the real and imaginary parts of the 
anomalous impurity Green function (G d ' ) are the same 
as those of the average anomalous Green functions of the 
superconductors. However, when the bound state moves 
below Ep, we find that the anomalous impurity green 
function changes sign. We interpret this to mean that 
the impurity has acquired an additional phase of ir. This 



behaviour is shown Fig. [|. It is clear from Eq. |T^ that 
this additional factor of (—1) will cause a sign change in 
the Josephson current. A similar phase shift wag found in 
earlier work involving a single superconductorQ, though 
a transition was not observed. 
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FIG. 4. Anomalous impurity spectral function ImG d ' 
for various values of the phase difference 8 with fixed 
<f>L + 4>r = 0, A = QMTk, T = 0.5 A. Note that as S in- 
creases, the spectral function changes sign; this coincides with 
the bound state crossing below the Fermi energy (see inset of 
Fig- 1). 

Conclusions.- Using an extension of the NCA, we have 
studied the properties of an infinite-!/ Anderson impu- 
rity coupled to two superconductors. We find that the 
sub-gap bound state develops an anomalous phase de- 
pendence, and that 7r-junction behaviour is lost when 
the bound state crosses above Ep. 
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